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In the paper we consider heuristic algorithm for solving graph isomorphism 
problem. The algorithm based on a successive splitting of the eigenvalues of 
the matrices which are modifications (to positive defined) of graphs' adjacency 
matrices. Modification of the algorithm allows to find a solution for Frobenius 
problem. Formulation of the Frobenius problem is following one. Given a pair 
of two matrices with the same number of rows and columns. We must find out 
whether one of the matrix can be acquired from another by permutation of 
it's rows and strings or not. For example, solution of Frobenius problem can 
give to us efficient way for decrypting of double permutation cyphers problem 
for high dimension matrices. 

The graph isomorphism problem 



Graph isomorphism problem is one of the problems for which we can't say def- 
initely whether this problem is polynomial or not [1]. But it's known that this 
problem is polynomial for some classes of graphs such as plane graphs, regular 
graphs and some others [2] , [3] , [4] . Our algorithm is heuristic algorithm for solving 
of the problem. 

We consider heuristic algorithm for solving the graph isomorphism problem. 
The algorithm based on successive splitting of eigenvalues of the matrices that 
are modifications (to positive defined) of graphs' adjacency matrices. During the 
work of the algorithm we solving the linear equations that defines inverse matrices 
of these matrices. Solutions of these systems gives to us a permutation which is 
sought bijection. 

There are two nonoriented graphs in the graph isomorphism problem: graph 
Ga — (VatEa) and graph Gb = {Vb,Eb)- Va, Vb are sets of graphs' vertices. 
Ea: Eb are sets of it's edges. They are supposed to be the sets of the same 
power: \Va\ = \Vb\, \Ea\ = \Eb\- Graph isomorphism problem formulation is 
following one: whether exists such bijection : Va Vb that if £ Ea then 
{<f{i),(f{j)) e Eb or not? 

Algorithm works with modifications of adjacency matrices. Let Aq be an adja- 
cency matrix of Ga , that is to say that — [o-lj ) and 

^0 ^ fl' if ihi) e Ea, 
[O, else. 

Let Bq be an adjacency matrix of graph Gb- 
Construct matrix Dao according to Aq: 

(di ... 0\ 
d2 ■•• 

VO ... dn) 
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Daq is a scalar matrix with following elements: 

n 

Construct matrix Dbq similarly according to Bq. Let 

A = Ao+Da„ B = Bo + Dbo. (1) 

Matrices A and B arc further considerated matrices. They are matrices which 
algorithm works with. They are symmetric and positive defined. 

If graps Ga = {Va, Ea) and Gb = {Vb, Eb) are isomorphic then we can acquire 
matrix A from matrix B by successive permutations of it's rows with simultane- 
ous permutations of it's columns with same numbers. So we can formulate the 
graph isomorphism problem as a particular case of the Frobenius problem which 
formulation is following one. Can we acquire matrix B from matrix A by successive 
permutations of it's rows and columns? 

Permutation of rows numbers i and number j for arbitrary matrix is equal to 
right multiplying matrix by permutation matrix P^ . Permutation of columns is 
equal to left multiplying matrix by the same permutation matrix. Permutation 
of vector's components with numbers i and j take place for both left and right 
multiplying vector by permutation matrix Pij . 

Let us consider the following systems of linear equations: 

Ax ^ ej,By = ek- (2) 

The vectors — (0, . . . , 0, 1, 0, . . . ,0) are basis vectors in the space i?", matrices 
A and B such as described above. Both systems are solvable and each solution is 
unique solution because of A and B are matrices with diagonal predominance and 
their determinants aren't equal to zero. Let be a solution of the system Ax = Cj 
and y'^ be a solution of the system By = Cfc. 

Note that solution of the systems of linear equations (1) gives to us inverse 
matrices of matrix A and matrix B. So for i-th component of vector x^ the following 
is true: x^ = (-l)*+^.4y7 det(A), where Aij are algebraic adjuncts for element (lij 
of the matrix A. That is to say x^ , j = 1, . . . ,n are columns of inverse matrix for 
A. 

If B = PjkAPjk then the following is true for solutions of the systems of linear 
equations (2): 

Xi = Vi, i^j, iy^k; 

Indeed: {Ax = ej) ^ {PjkAx = Pjkej) ^ {PjkAxPjk = PjkejPjk) ~ {PjkAPjkX 
■ ~ {PjkAPjkxPjk = ejPjk) ~ {BxPjk = efe). That is to say xPjk = y- 

If matrix B can be acquired from matrix A by successive permutations of it's 
rows with simultaneous permutations of it's columns with same numbers then: 

B = Pjiki ■ ■ ■ Pjiki-^Pjiki • • • Pjikii 

and consequently xPj^ki ■ ■ ■ Pjiki = V- 

So if we'll change vector Cfc in the system of linear equations (2) with fixed j 
by changing index k from 1 to n then vectors x^ and y^ ^ which are corresponding 
solutions of the systems (2), will be the same vectors correct to permutation of 
their components only if row with number j of matrix A corresponds to row with 
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number k of matrix B. That is to say elements of row number k of matrix B are 
permutatcd elements of row number i of matrix A. The same is true for columns 
of the matrices. 

We'll accomplish the successive perturbations of its diagonal elements for finding 
unique row and unique column of the matrix B that is corresponds to row and col- 
umn with number i of the matrix A. The algorithm works with symmetric matrices 
which can be transformed into scalar matrices by orthogonal transformations. So 



where A, B are scalar matrices with eigenvalues on their diagonals and Ua, Ub 
are matrices of the orthogonal transformations. Diagonal elements of A and B are 
eigenvalues of the matrices A and B. 

The spectrums of the isomorphic graphs' adjacency matrices are the same [5]. 
The same true for the matrices which algorithm works with because of described 
above procedure of changing of the diagonal elements of adjacency matrices leads 
only to shifting spectrums of the matrices. So if the matrices are corresponds to the 
isomorphic graphs then spectrums of the matrices will be congruent after shifting. 

It's obvious that if the spectrums of the both matrices is simple or multiplicity 
of eigenvalues is not big then the problem of graph isomorphism can be solved by 
a comparison of rows of the matrices J7a, Ub, A and B [7]. 

The main difficulties arise with consideration of graphs which spectrums contains 
multiple eigenvalues. Perturbation of the matrices will allow us to perturb their 
spectrums while algorithm works. Splitting of the eigenvalues will take place by 
this perturbation. So we can maintain an unique correspondence between rows and 
columns of the matrices. 

If there is multiple eigenvalues in spectrums of the matrices and then they will 
be split and we'll be able to find the permutation which is sought bijcction Lp. This 
bijection establishes isomorphism of graphs Ga — (Vk, Ea) and Gb — (Vs, Eb)- 
Computing experiments give us evidence that split needed for determination of the 
corresponding rows and columns take place much early than at the last iteration 
of the algorithm. So at iteration with number ^Jn there is no need for further 
perturbations of diagonal elements of the matrices A and B for maintaining a 
unique correspondence between vertices of graphs that represents lattice on torus 
which numbers of vertices is equal to the n where n changes from 9 up to 400 . 

At the algorithm's implementation on every iteration algorithm works with al- 
ready perturb matrices not with initial matrices. So if we've got a correspondence 
between row number j of matrix A^ and row number k of matrix B^ at the iteration 
number j and for columns with the same numbers then we'll considerate further 
perturb matrices A^'^^ and B^~^^: 



We make perturbation by scalar matrices C'^ with elements Cj on it's diagonal: 



As a result if matrix A can be acquired from matrix B by successive permutations 
of it's rows with simultaneous permutations of it's columns with same numbers then 



A = UaAUI, 
B = UbBUI, 



A^+^ = A^ +eC^, B^+ 



B^ +eC''. 




1, if i = j = k 
0, else. 
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we get the sought permutation P while j changes from 1 to n. That is to say 

P=(^ 2 ... n 
\ki k2 ... fc„y 

kj is a number of row of matrix B tliat obtained at j-th iteration of the algorithm. 

As a matter of fact permutation P is one of the possible permutations. So P is 
a bijection <fi :Va ^ Vb that sets isomorphism of graphs Ga and Gb- 

The spectral splitting algorithm 

(First schema) 

Step 0. ^0 := A, j := 1. 

Step 1. If j <n then := A^~^ + eC^ , else stop the algorithm's implementation. 
Step 2. Solving of the system of linear equations Ax = ej. xJ is the solution. 
Step 3. fc := 1. If < n then go to Step 3.1, else go to Step 4. 
Step 3.1. B'' :=B'=-i +£C'=. 

Step 3.2. Solving of the system of hnear equations B''y = Cfc. y'^ is the solution 

Step 3.3. k:=k+l. Go to Step 3. 

Step 4. Comparing norms of and y'^ , k = 1, . . . ,n. 

If Vfc I jx-' 1 1 ^ I Ij/*^ 1 1 then graphs Ga and Gb are not isomorphic. Stop the algorithm's 
implementation. 

If exists A; : ||a;-'|| = = then P{j) := k (Maintaining a correspondence 

between vertex j of graph Ga and vortex k of graph Gb), := B^^^ + eC^ . Else 
graphs Ga and Gb are not isomorphic. Stop the algorithm implementation. 
Step 5. j := j + 1. Go to Step 1. 

The hardness of this scheme is equal to O(n^), n is a number of rows at the 
square matrices A and B. We have to notice that this scheme can be modified to 
another scheme which laboriousness is equal to 0{ii?-^). We can make it by adding 
the split checking procedure (Step 6 and step 7 of the scheme stated below). 
The spectral splitting algorithm 
(Second schema) 

Step 0. AO := A, j := 1. 

Step 1. If j < n then A^ :— A^~^ + eC^ , else stop the algorithm's implementation. 
Step 2. Solving of the system of the linear equations Ax = ej. x^ is the solution. 
Step 3. k := I. If fc < n, then go to Step 3.1, else then go to Step 4. 
Step 3.1. B'^ BJ-i +£C'=. 

Step 3.2. Solving of the system of the linear equations B'^y = Ck- y'^ is the solution. 

Step 3.3. k:= k + 1. Go to Stop 3. 

Step 4. Comparing norms of x^ and y'' , k — 1, . . . ,n. 

If Vfc I Ix-' 1 1 7^ I ly'^ 1 1 then graphs Ga and Gb are not isomorphic. Stop the algorithm's 

implementation. 

If exists k : \\x^\\ = \ \y'^\\ Xj = y^ then P{j) := k (Maintaining a correspondence 
between vertex j of graph Ga and vertex k of graph Gb), := B^~^ + sC'^. Else 
graphs Ga and Gb are not isomorphic. Stop the algorithm's implementation. 
Step 5. j := j + 1. 

Step 6. Solving of the system of linear equations A^x = ei. x', Z = 1, . . . , n is the 

solution. 

Step 7. If VZ yp:\\x^\\^ \ \xP\\ Go to Step 8, else go to Step 1. 

Step 8. If j < n then implement solving of the system of the linear equations 
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B'^y = Cfc, fc = 1, . . . , n; where k such that not exists i : P{i) = j. 

Step 9. Comparing norms of and Wy'^W, I = j, . . . ,n, k = 1, . . . ,n; and k such 

that not exists i : P{i) =j. If = ||/|| then P{1) := k. 

Step 10. Stop the algorithm's implementation. 

Though the spectral splitting algorithm show oneself good by computing ex- 
periments we have to notice that we can't guarantee finding solution of the graph 
isomorphism problem for arbitrary graphs because of we can't say that the follow- 
ing situation is impossible. There are two vectors Xj and jjk with equal or such 
close norms for we can't say whether one vector can be acquired from another by 
permutation of it's components or not because of we implements check procedure 
only by one components of each vectors. They are components Xj and y^. 

In particular we realized computing experiments for regular graphs with number 
of vertices up to 2500 for the graph isomorphism problem. 



Solving of the Frobenius problem for Eirbitrary 
squcire matrices of complete rank 

The spectral splitting algorithm for the graph isomorphism problem can be ap- 
plied for weighted graphs without any modifications. Scheme of the algorithm is 
the same; the only difference in matrices algorithm works with. 

There are two nonoriented graphs in the graph isomorphism problem: graph 
Ga = {Va,Ea) and graph Gb = {Vb,Eb)- Va, Vb are sets of graphs' vertices, 
Ea, Eb are sets of it's edges. They are supposed to be the sets of the same power: 
\Va\ = \Vb\, \Ea\ = \Eb\. Given a function Ha : Ea ^ R Hb : Eb ^ R 
which defines weights of graphs' edges. Graph isomorphism problem for weighted 
graphs formulation is following one: whether exists such bijection ip : Va ^ Vb, 
that if e Ea, then {ip(i),ip{j)) e Eb and if e Ea then , ip{j)) e Eb 
HA{i,j) = HB{ip{i),ip{j)) or not? 

The adjacency matrices of weighted nonoriented graphs transforms into positive 
defined matrices with diagonal predominance, but we should select its diagonal 
elements in such way that conditionality number should be confined. 

Let matrix Aq = (a°j) be an adjacency matrix of weighted nonoriented graph 
Ga and matrix Daq be a diagonal matrix with such elements di that: 

n 

di = d + ^a°j, 

n 

where d = max ^ a? ■ . 

l<i<7i j—l 

Construct matrix Dbq for weighted nonoriented graph Gb by same way accord- 
ing to adjacency matrix of graph Bq. 

For conditionality number IJl{A) of the symmetric matrix A the following is true 
[6]: 



where ri{A) = max {an -\- X] \aij\), x{A) = niin {an - X] kijl)- 

1<»<" j^i l<i<n 

By our choice of d: 

7~t 71 71 Tt 

r,{A) = max {au + E K\) = «nn + E Kj\ = c^n + E Kj\ = ^ + 2 E Kj\ 

!<»<" , = 1 .,=1 ,-=1 ,-=1 
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i=i 

n 

X{A) = mm {an - E |ai, l) 



l<i< 



i=i 



E l< 



El<,l + rf- 



- E K,\ = E l< 

Consequently 



■EI<,|-EI<,I = EI<,I- 
j=i j=i j=i 



= 3. 



There are two matrices Fa and Fb with the same number of rows and columns 
in the Frobenius problem. We must find out whether one matrix can be acquired 
from another by some permutation of it rows and columns. 

Suppose that matrices Fa and Fb are square matrices of the complete rank with 
number of rows is equal to n. Matrix of the following structure: 

Fa'' 
FJ 

This matrix can be considered as an adjacency matrix of some weighted nonoriented 
graph. Construct matrix Daq with structure described above for matrix Aq. Then 
construct matrix A: 

A = Ao + Dao. 

That is to say 



Fa 



and 



where 



fdi 






d2 



\0 



Dao 



0\ 





dnj 



( dn+1 





V 





dn+2 





\ 





d2n/ 



So matrix A is matrix for initial matrix Fa for algorithm to work with. 

Matrix A is positive defined symmetric matrix with diagonal predominance. It 
corresponds to some weighted nonoriented dicotyledonous graph. Matrix Fb is a 
square matrix which number of rows the same as the number of rows of matrix Fa- 

Construct matrix B by the same way for matrix Fb- Form the solving permu- 
tation of rows and columns of matrices A and B by applying the spectral splitting 
algorithm to matrices A and B- 

Form the solving permutation for initial matrices Fa and Fb by the permutations 
of rows and columns of the matrices A and B in the following way. Let it given by 
the algorithm that the row and column with number i of matrix A corresponds to 
the row and column with number j of matrix B. It is obvious that if 1 < i < n, 1 < 
j < n then the permutation of pair of rows with simultaneous permutation of 
it's columns {i,j} of matrix A corresponds to permutation of rows number i and 
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number j of initial matrix Fa- Ifn+l<i<2n, n+1 < j <2n then permutation 
of pair of rows {i,j} with simuhancous permutation of columns {i,j} of matrix 
A corresponds to permutation of columns with number i and number j of initial 
matrix Fa- As soon as matrices A and B which algorithm works with have the 
structure as following one: 



/* 
* 





* * 



* * 
* * 



V* 



* * * 

* * 

* * 

* : : 

* 



(* are positions admissible for nonzero elements) so the situation of maintaining of 
the correspondence for pair of rows and pair of columns of matrices A and B with 
such numbers i and j that {i, j} and 1 < i < n, n + 1 < j < 2n is impossible. It 
impossible because of it means appearance of nonzero nondiagonal elements at the 
left and upper part of the matrix A (of the matrix D\^) and right and lower part 
(of the matrix D'j^^) if only the initial matrices Fa and Fb are not scalar. This is 
impossible by constructing of the matrices A and B. Consequently appearance of 
such correspondence is impossible while algorithm works. 



Deciphering of the double permutation cipher 
by the spectral splitting algorithm 

Suppose we have a text presented by some square matrix of symbols. Let us 
form matrix of the cipher with some digital code for every symbol of text. The 
second matrix presented some coding by this way text too. It's needed to find out 
whether the text presented by the second matrix is the double permutation cipher 
for the firs matrix or not. That is to say whether this text can be obtained from 
the first text presented by matrix by permutation of it's rows and columns. Such 
problem may be interpreted as a case of the Frobenius problem for digital matrices. 
So we can apply the spectral splitting algorithm for solving of this problem. 

There are computational experiments was implemented for deciphering of the 
double permutation ciphers for texts with number of symbols up to 10000. 

References 

1. M. Garey and D. Jhonson Computers and Intractability, A Guide to Theory 
of NP-completness. Freeman and Co., 1979 

2. Hopcroft, Wong A linear time algorithm for isomorphism of planar graphs. 
Proceedings of the Sixth Annual ACM Symposium on Theory of Computing, p. 
172-184, 1974. 

3. Luks Isomorphism of graphs of hounded valence can he tested in polynomial 
time. Proc. 21st IEEE FOCS Symp., 42, 49, 1980. 



8 



4. Hoffmann Group- Theoretic Algorithms and Graph Isomorphism. Lecture 
Notes in Computer Science (Chapter V). P.127-138, 1982 

5. Tzvetkovich D. and others Spectrums of graphs. Theory and application. 
Kiev: The scientific thought, 1984. 

6. Godunov S. and others Guaranteed exactness of solutions of systems of linear 
equations in euclidian spaces. Novosibirsk: Science, 1988. 

7. Kikina A., FaizuUin R. The algorithm for testing graph isomorphism. VINITI 
21.06.95 1789-95. 



